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ABSTRACT 

We present an analysis of the neutral hydrogen and stellar populations of elliptical galaxies in 



the Tal et al. (2009) sample. Our aim is to test their conclusion that the continuing assembly of 



these galaxies at z~0 is essentially gas-free and not accompanied by significant star formation. 
In order to do so, we make use of H i data and line-strength indices available in the literature. 
We look for direct and indirect evidence of the presence of cold gas during the recent assembly 
of these objects and analyse its relation to galaxy morphological fine structure. 

We find that >25 % of ellipticals contain H i at the level of M(H i)> 1 8 M G , and that M(H i) 
is of the order of a few percent of the total stellar mass. Available data are insufficient to es- 
tablish whether galaxies with a disturbed stellar morphology are more likely to contain H i. 
However, H i interferometry reveals very disturbed gas morphology/kinematics in all but one 
of the detected systems, confirming the continuing assembly of many ellipticals but also show- 
ing that this is not necessarily gas-free. We also find that all very disturbed ellipticals have a 
single-stellar-population-equivalent age <4 Gyr. We interpret this as evidence that ~0.5-5% of 
their stellar mass is contained in a young population formed during the past ~1 Gyr. Overall, a 
large fraction of ellipticals seem to have continued their assembly over the past few Gyr in the 
presence of a mass of cold gas of the order of 10% of the galaxy stellar mass. This material is 
now observable as neutral hydrogen and young stars. 
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1 INTRODUCTION 

The role of dissipationless galaxy merging in the assembly of early- 
type galaxies is currently under debate. On the one hand, it ap- 
pears to be necessary to produce the boxy, slowly-rotating galaxies 
which populate the high-mass end of the red sequence (e.g., |Faber| 
et al. 1997; Kormendy et al. 2009 and references therein). On the 
other hand, dissipative processes are required to explain a num- 
ber of structural, kinematical and stellar-population properties of 
intermediate- and low-mass early-type galaxies (e.g., [Bender et al. 
[19921 |Emsellem et"aT1|2007| |Hopkins et al][2009| and refernces 
therein). Furthermore, the tight scaling relations observed in the lo- 
cal Universe place a rather conservative upper limit on the fraction 
of stellar mass assembled via dissipationless merging (e.g., |Nipoti| 
et al. 2003, 20091. Finally, cold gas is actually observed in a large 
fraction of early-type galaxies in the nearby Universe (e.g., |Mor-| 
|ganti et al.|2006||Combes et al.|2007] >. 

Based on deep imaging of a volume-limited sample of nearby 
elliptical galaxies with M s ^ -20, and following earlier work by, 
e.g.,|Malin & Carter| <|T983]> and |Schweizer & Seitzer| < [T992] >, [Tall 



|et al. ( 2009 hereafter T09) find signatures of recent dynamically 
violent assembly in a remarkably large fraction of objects, ~75%. 
Analysing galaxies' B-V colour, they argue that ellipticals outside 
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cluster environments continue to grow at z=0 mostly through gas- 
free accretions. If correct, this conclusion would have implications 
for a number of properties of local ellipticals. This follows from 
a large number of numerical studies showing that the presence of 
even a modest amount of gas (~ 10-20% of the baryonic mass) in 
a binary merger can change the stellar mass-distribution and orbits 
in the early-type remnant relative to the gas-free case (e.g., |Jes-| 
|seit et al.|2005l|2007||Hopkins et al.|2009| >. Observable quantities 
such as ellipticity, isophote disciness/boxiness, anisotropy and light 
profile are affected by the gas content of the progenitor galaxies. 

Although most of these numerical simulations focus on binary 
major or minor mergers, early-type galaxies in a ACDM Universe 
form following a sequence of many interactions with widely vary- 
ing mass ratio and degree of dissipation. In other words, galaxies 
are much more likely to grow via continuous assembly of smaller 
units of varying gas content rather than a few major events, in a 
way strongly dependent on their environment (e.g., |De Lucia et al.| 
[20561 |Naab e t al.|2007||Parry et al.|2009) . Within this framework, 
as in the binary merger case, a galaxy formed following a gas-free 
path (e.g., one residing in the centre of a massive cluster, where 
little cold gas can penetrate) should be significantly different from 
one that experienced dissipation (e.g., Bournaud et al. 2007). It is 
therefore interesting to estimate the mass of cold gas present during 
the continuing assembly of ellipticals over the last few Gyr and test 
T09 conclusion that this is indeed negligible. 

The approach taken in this Letter is to look for neutral hydro- 
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Figure 1. H i content of galaxies in the T09 sample. Circles are H i-detected 
galaxies, triangles are upper limits on the Hi mass. Colour represents the 
environment as reported by T09. 



gen and young stellar populations in galaxies in the T09 sample. 
Both properties could be regarded as marginally important with re- 
spect to the broad cosmological picture of galaxy evolution. After 
all, work over the past decades shows that young stars and cold gas 
amount to just a few percent of the total galaxy stellar mass of z=0 
ellipticals. However, as explained above, even modest amounts of 
gas can have important effects on fundamental properties of these 
objects. We would like to clarify what kind of data analysis is nec- 
essary in order to detect this gas. Work along this line has been 
carried out by a number of authors during the past few years. We 
summarise it together with our results in the discussion section. 



2 NEUTRAL HYDROGEN 

All galaxies in the T09 sample have been observed in the 21 -cm 
line by the HIPASS survey ( |Meyer et al.|2004) . The noise in the 
HIPASS spectra is ~15 mjy in ~13-km/s-wide channels. It cor- 
responds to a 3cr flux of ~2.4 Jy-km/s over 200 km/s, which we 
take as the typical width of the H i profile. Therefore, the HIPASS 
M(Hi) detection limit varies in the range 0.1-1.2xl0 9 M Q as a func- 
tion of distance for the T09 sampl^] This is not very deep and in- 
deed, after inspection of the HIPASS spectra, we find that only 2 
galaxies are detected: IC 1459 and IC 4889. We searched the liter- 
ature for deeper Hi observations. This resulted in 12 more detec- 
tions and 9 upper limits deeper than those provided by HIPASS 
(derived with the same criterion described above). In one more 
galaxy, NGC 6868, H i is detected in absorption against the central 
radio continuum source (Oosterloo, priv. comm.) and estimating its 
M(H i) is not possible with current data. We refer to Table [TJ for 
details on the H i data. 

In Fig[TJ we plot M(H i) and M(H i)/L B against the tidal pa- 
rameter T c introduced by T09. T c is designed as a way of quanti- 
fying the amount of morphological fine structure and is larger for 
galaxies with higher disturbance. In the figure, colour represents 
galaxy environment following the T09 division of the sample in 



cluster, group, field and unknown environment. Overall, we find 
Hi in/around ~l/4 of ellipticals, and mostly outside the cluster en- 
vironment. Because of the inhomogeneity of the data, we conclude 
that >25% of ellipticals contain Hi at the 10 8 -M Q level. M(Hi) can 
be as high as a few times 10 9 M Q , corresponding to at most a few 
percent of the total stellar mass. Given the many upper limits, we 
cannot establish whether highly-disturbed ellipticals host more H i 
in/around their stellar body than more relaxed objects. 

Most of the H i data used here are taken with single-dish tele- 
scopes. These have a beam of several arcmin FWHM and we have 
discarded 4 detections that are likely caused by the presence of a 
late-type neighbour within the same beam. However, a number of 
observations taken with interferometers are available and allow us 
to study the actual spatial distribution of the H i around about half 
of the detected galaxies. Here we comment on individual galaxies 
for which radio interferometry is available. 

IC 1459 shows a number of Hi filaments surrounding the stel- 
lar body; the galaxy resides in a rich group with many gas-rich disc 
galaxies. IC 4889 hosts a strongly-warped disc/ring, but the poor 
spatial resolution of the ATCA data may hide a more complicated 
configuration. NGC 2865 is a known recent-merger remnant with 
H i associated to a number of stellar shells. NGC 2974 hosts a reg- 
ular H i ring aligned with the stellar body. NGC 4472, the bright- 
est galaxy in the Virgo cluster, shows an Hi cloud that suggests 
on-going interaction with a dwarf neighbour. NGC 5018 shows a 
filament of gas stretching across its stellar body and connecting 
two neighbours. NGC 5077 is likely in the process of acquiring gas 
from a late-type companion as a broad tail of H i connects the two 
systems and H i clouds are scattered around them. NGC 5903 hosts 
kinematically-disturbed H i distributed along the galaxy minor axis 
and two gas tails connected to the edge of this distribution. 

Overall, interferometric data show a disturbed H i configura- 
tion, signature of on-ongoing or recent (~1 Gyr) gas accretion, in 
all but one of the detected ellipticals. Considering gas consumption 
in star formation and tidal-tail gas lost to the inter-galactic medium, 
it is plausible that many ellipticals have continued their assembly 
until recently in the presence of a mass of cold gas as high as ~10% 
of their total stellar mass. Given the many high upper limits on 
M(H i), deeper H i observations would be required to establish the 
actual fraction of ellipticals that host neutral hydrogen at this level 
in the T09 sample. However, we can state that this is >25%, in 
agreement with deep observations where H i is found in about half 
of the targets (Morganti et al. 20061. It is clear that the presence of 
significant amounts of cold gas during the recent assembly of these 
systems is fairly common outside cluster environments. 



3 STELLAR POPULATIONS 

Fully corrected line-strength indices measured from optical spectra 
and brought onto the Lick/IDS system are available in the literature 
for 38/55 galaxies in the T09 sample. For 34 of them indices within 
an RJ10 aperture are given by Thomas et al. (20051 or | AnnibalT| 
|et al.| ( |2006^ . Indices for a few additional objects and measured 
over slightly different apertures are available in Sanchez-Bl azquez] 
[eTaT1p007t and |Serra etaT1p008} . We refer to Table[TJfor details 
on the line-strength data. 

In Fig[2]we plot H/?, an age-sensitive index, against [MgFe], a 
metallicity indicator nearly insensitive to [ff/Fe] variation^ Galax- 



M(Hi)=2.36 X 10 5 F(Hi)/(Jy ■ km/s) (d/Mpc) 2 M e 



[MgFe]= ^Mgb X (Fe5270 + Fe5335)/2 
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Figure 2. H/3 vs. [MgFe] for the 38 galaxies in the T09 sample with avail- 
able line-strength indices. Colour codes the tidal parameter T c defined in 
T09. The model grid represents SSP models from |Thomas et aT|j2003f . 
Solid lines are models of fixed metallicity with [Z/H] labelled at the bot- 
tom of the grid. Dashed lines are models of fixed stellar age with. 
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Figure 3. fssp vs. r c for the 38 galaxies in the T09 sample with available 
line-strength indices. The grey-scale codes Mg. 



ies are plotted on top of a model grid which represents solar-fa/Fe] 
single-stellar-population (SSP) models of Tho mas et al.| (20031. 
Colour represents T c . The figure shows that high-T c ellipticals tend 
to be younger than more relaxed objects. In particular, there are no 
old galaxies with r c >0.13. The galaxy with very large central HB 
is NGC 2865, a known gas-rich merger remnant. 

We derive best-fitting SSP-equivalent age and metallicity (fssp 
and [Z/H] SSP ) by comparing line-strength indices to predictions 
from |Thomas et al.] ( |2003[ > models. The comparison is performed on 
the [MgFe]-H/3 plane shown in Fig[2] We have verified that fssp, the 
most relevant parameter for our study, is in agreement with values 
taken from Thomas et al. (20051, Annibali et al. (20071, Sanchez- 
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(20081 despite the fact that 



Blazq uez et al.| ( [2"007[ l and 

we neglect the effect of [a/Fe]. We list fssp values in Table [T] 

In Fig[3]we plot f S sp against T c , with grey-scale representing 
M B . While we do not see a clean correlation between fssp and T c , 
highly-distrubed ellipticals clearly stand out as they are all fairly 
young (f SS p<4 Gyr). On the contrary, most galaxies with low T c 



are older, although there are a number of weakly-disturbed objects 
with young ages too. We note that high-r c galaxies are not younger 
because of being systematically fainter than more relaxed objects. 

If a galaxy hosts stars with a spread in age, fssp is strongly bi- 
ased towards the age of the youngest stars (see |Serra & Trager|2007| 
for a thorough discussion of this effect). Therefore, it is likely that 
galaxies with a fssp of just a few Gyr host recently-formed stars on 
top of an old, dominant population. We estimate the young-to-total 
stellar mass fraction M y /M t by comparing HB and [MgFe] to pre- 
dictions from solar- [a/Fe], two-SSP models built based on Thomas 
et al. (20031 SSPs. Since the number of free parameters is large, 
we assume that the old population has an age of 13 Gyr and that 
the young population has solar metallicity. We let the metallicity of 
the old population (which dominates [MgFe] - see |Serra & Tragerj 
|2007[ > and M y /M t free to vary, and fix the age of the young pop- 
ulation f y to 300 Myr and 1 Gyr in two different fitting iterations. 
Given the known degeneracy between f y and M y /M t this approach 
allows us to more easily estimate the allowed range for M y /M t . 

We find that ellipticals with fssp<4 Gyr contain 0.5-5% of their 
stellar mass in stars formed 300 Myr to ~ 1 Gyr ago. Given the many 
approximations involved in this analysis these numbers should be 
taken with some care. Nevertheless, it is hard to question that many 
ellipticals, and in particular all highly-disturbed objects, are con- 
sistent with having formed up to a few percent of their stellar mass 
within the past ~ 1 Gyr. We note that it would have been impossible 
to reach this conclusion based on optical photometry only. Indeed, 
the best-fitting two-SSP models predict B - V colours in the range 
0.9-1.1 with typical lcr uncertainty below 0.1 dex. The only excep- 
tion is NGC 2865, which is predicted to have B - V -0.8. 



4 DISCUSSION AND CONCLUSIONS 

Neutral hydrogen observations and line-strength index analysis of 
the stellar populations of ellipticals in the T09 sample reveal that a 
large fraction of these galaxies continued their assembly during the 
past few Gyr in the presence of a mass of cold gas of at least a few 
percent of the total galaxy stellar mass. This gas is now detected in 
the form of Hi and through young stellar sub-populations. 

H i interferometry reveals that in nearly all cases the detected 
gas is morphologically/kinematically disturbed and mostly found in 
tails/filaments, sometimes clearly associated to the optical morpho- 
logical disturbances. This, and the presence of young stars, confirm 
T09 conclusion that the assembly of elliptical galaxies continues to 
z=0. However, it does not agree with their conclusion that this as- 
sembly is essentially gas-free. In fact, Hi observations could be 
regarded as yet one more tool to reveal the on-going assembly of 
these systems, complementary to deep optical imaging (see |Mor^| 
|ganti et al.|2006)|Oosterloo et al.|2007l|Grossi et al.|2009> . 

Our estimate of the mass of cold gas recently accreted by el- 
lipticals is in qualitative agreement with the mass accretion rate 
derived by T09 themselves. They estimate that nearby ellipticals 
are still growing at z=0 at a rate of 20%-in-mass per Gyr. Evidence 
of this on-going assembly is particularly strong in the field and in 
galaxy groups, where most small galaxies around a given elliptical 
are gas rich. It is therefore reasonable to expect that such a conspic- 
uous mass assembly is accompanied by an accretion of cold gas of 
at least a few percent of the total galaxy stellar mass. 

Despite this high incidence of disturbed Hi systems, ellip- 
tical galaxies with very regular gas distributions do exist (e.g., 
NGC 2974; more examples can be found in Morgant i et al.|2006} 
|Oosterloo~et al. 2007 1. These systems could be evolved versions of 
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the galaxies that show disturbed Hi at z=0, i.e., galaxies for which 
gas accretion was important many Gyr ago and whose recent his- 
tory has been quiet enough for this gas to settle. 

Another interesting point is that while we find that optically 
disturbed ellipticals are systematically younger than relaxed ob- 
jects, and interpret this as evidence of their continuing, gas-rich 
assembly, not all ellipticals with Hi are young (e.g., NGC 5846; 
other cases are known outside the T09 sample, e.g., NGC 4278 in 
Morg anti et al.|20 06 ). Clearly, the dynamics of the accretion plays 
a very important role in determining the effect of gas on structure 
and (distribution of) stellar population of the accreting galaxy. Un- 
derstanding where the young stars and cold gas are, and what their 
kinematics is, is as important as estimating their mass. 

Our results agree with previous work on this subject. [Sanchez-| 
Blazquez et al. ( 2009 ) analyse the stellar population of high-z early- 
type galaxies selected by |van Do kkum (2005 ) to be "dry"-merger 
remnants. They find that highly disturbed objects are on average 
younger than relaxed ones, and that their line-strength indices are 
consistent with the presence of a few-percent-in-mass of young 
stars. This is in agreement with early results from Schweizer et al. 
( 1990) who find the Hy3 index to be systematically higher for el- 
lipticals with morphological fine structure at any given absolute 
magnitude. Finally, Donovan et al. ( 2007 1 look at the H i properties 
of a sample of local galaxies selected in the same way as the [van] 
Dokkum ( 2005 1 objects. They find that many of these systems host 
a significant mass of H i whose disturbed morphology and kinemat- 
ics make them consistent with having recently accreted cold gas. 

Following the conclusions of a large number of numerical 
studies, the direct or indirect detection of gas present during the 
assembly of elliptical galaxies may have important implications for 
many structural properties of these objects. We have shown that us- 
ing only observables such as colours measured over large apertures 
one risks to neglect the evidence of this gas because of the lower 
sensitivity to the presence of a modest, but nevertheless significant, 
mass of cold gas or young-stars. 
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(1) Distance from: (a) Tonry et al. (2001 1; (b) NED recessional velocity corrected for Virgo infall (h = 0.73); (c) NED z-independent. (2) Mb derived from 
HyperLeda btc. (3) Colour as in T09. (4) T c from T09. (5) M(Hi) adopting the total Hi flux from: (d) Oosterloo (priv. comm.); (e) HIPASS spectra (Meyer 
et al.|2004F; (/ ) |Oosterloo et al.|(2007); ( g) |Haynes^TaT]{l99 7 j); (/i)|Sansom et al.H2000>; (Q|Knapp et aTHl985); (;)|Huchtmeier|(l994); (^jSchiminovich 
et al.|l '1995 1; (/) Weijm ans et al.H2008) ; (m) Bottinelli & Gougue nheimH1977); (n)|Bottinelli & G ouguenhe im|(1979b); (o)|Bottinelli & Go uguenheim 
(1979a> ; (p) |McNamara et al.|(199lT rto) |Knapp et al.H1978) ; (r) |Kim et al-H1988) ; (^) |Lake & SchommerH1984t ; (Q |LewisH1987t ; (M) |Appleton et"aL| 
41990); s ee text for the calculation of the upper limits. (6-8) Line-strength indices from: (v) ]ATm1)a!i"eta!*|j2006F ; (w) Trager (priv. comm.); (;t) |Thomas et al. 
(2005) ; Cy) |Serra et alj(2008) ; (z) |Sanchez-Blazquez et aTH2007) ; <Fe>=(Fe5270+Fe5335)/2. (9) fssp derived in this work. (10) Environment from T09. 
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